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Abstract 

The isoscaling parameter usually denoted by a depends upon both the symmetry energy coeffi- 
cient and the isotopic contents of the dissociating systems. We compute a in theoretical models: 
first in a simple mean field model and then in thermodynamic models using both grand canonical 
and canonical ensembles. For finite systems the canonical ensemble is much more appropriate. 
The model values of a are compared with a much used standard formula. Next we turn to cases 
where in experiments, there are significant deviations from isoscaling. We show that in such cases, 
although the grand canonical model fails, the canonical model is capable of explaining the data. 

PACS numbers: 25.70Mn, 25.70Pq 
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I. INTRODUCTION 



For central collisions of Sn on Sn ( 112 Sn+ 112 Sn, 124 Sn+ 112 Sn and 124 Sn+ 124 Sn) a well- 
known result is that the ratio of isotope yields from two different reactions, 1 and 2, 
R 2 i(N,Z) = Y 2 (N, Z)/Yi(N, Z) exhibits an exponential relationship as a function of the 
isotope neutron number N and proton number Z 

R 21 (N, Z) = Y 2 (N, Z)/Y 1 (N, Z) = C ex V (aN + pZ) (1) 

This is called isoscaling. Note that for Sn on Sn central collisions the fragmenting system 
is rather large. When the fragmenting system is significantly smaller, the above equation is 
only approximate |4|. We will confine ourselves to large systems till we come to sections VII 
and VIII. 

Much effort has gone into trying to relate a to the symmetry energy term that occurs in 
liquid drop binding energy formula. In its simplest version the symmetry energy term is given 
by C S (N — Z) 2 /A. It is reasonable to guess that the ratio R 2 i(N, Z) should predominantly 
depend on N ,Z or equivalently on Z ,A Q of the fragmenting systems and also on the 
value of C s . An approximate functional relationship that can be deduced from models is: 



a 



4^((Z (1)A4 (1)) 2 - (Z (2)Mo(2)) 2 ) where Z (l)/A (l) refers to the disassociating 
system in reaction 1, Zq(2)/Ao(2) refers to that in reaction 2 and T is the characteristic 
temperature in the two reactions. It is this approximate equality that we examine in this 
work, first in a mean-field model (section III) and then, in detail, in thermodynamic model 
using both canonical and grand canonical ensembles. In particular we point out a different 
functional relationship is more natural in certain physical situations. Next we turn to cases 
where isoscaling is only approximate and how such cases can be handled in the theoretical 
framework. Summary and conclusions are presented in section X. 



II. RELATING ISOSCALING PARAMETER TO CHEMICAL POTENTIAL 

Equation (1) can be easily understood using a grand canonical model for multifragmen- 
tation. This allows us to relate a to chemical potentials. 

We assume that in a central collision, the two ions fuse, some pre-equilibrium emission 
occurs and the fused system, because of two-body collisions is heated up and begins to 
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expand. During the expansion composites are formed. As the expansion takes place inter- 
action between compsites rapidly fall off except for Coulomb interaction which can be taken 
care of in an approximate way using the Wigner-Seitz apporoximation In this expanded 
volume the break up of the dissociating system can be calculated using laws of equilibrium 
statistical mechanics. The calculation is particularly simple if a grand canonical ensemble 
is used. 

More about the grand canonical approximation will follow later but at this stage let us 
quickly connect a to chemical potential encountered in the grand canonical ensemble. The 
cross-section of the produced composite is given by a(N, Z) = C{riN^z) where C is a constant 
not provided by the model; (riN,z) is average multiplicity of the composite. For system 1 
characterized by total charge Z (l) and total mass Aq(1) (total neutron number N (l) = 
Aq(1) — Z (l)) and neutron and proton chemical potentials /x n (l) and // p (l) respectively this 
multiplicity is 

(n N>z (l)) = efi^m+fi^0)z UNtZ (2 ) 

where (3 is the inverse of temperature T and ujn,z is the one particle partition function of 
the composite N, Z. 

If in the second reaction the total charge is Zq(2), the total mass is Aq(2) but the con- 
ditions of the second reaction are similar to that of reaction 1 and we expect the same 
temperature, then 

as the uj(N, Z) in the numerator cancels the u(N, Z) in the denominator. Let Sfi n = /i n (2) — 
H n (l) and 5fi p = /i p (2) — A widely used relationship is 

0fin Cs[[ A (lY [ A (2) )l [) 
A corresponding relationship can be written down for 5/j, p . It suffices to study any one and 
traditionally one examines 5{i n . In the following we will investigate 5fi n for various cases. 



III. 6 ti n IN A MEAN FIELD MODEL 



The concept of chemical potential is useful not only in problems concerned with mul- 
tifragmentation. We first investigate the chemical potential in mean field theory at finite 
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temperature. One might argue that this is a valid model at low temperatures T < 3 MeV. 

n q 

The caloric curve has been computed in this model [g, |7| and many interesting results were 
found. 

For a nucleus with iV neutrons and Z protons (N + Z = A) the symmetry energy 
contributes to the binding energy a term: C s ■ To binding energy per particle it gives 

C S ^ N ^P ■ The term C s has its origin to both kinetic and potential energy per particle so 
we separate C s into two parts: C s = C s (k.e) + C s {p.e). 

We consider asymmetric nuclear matter where proton charges are switched off. We can 
calculate both C s (k.e) and C s (p.e) in the Hartree-Fock model interacting by Skyrme inter- 
action. For an asymmetric nucleus p n and p p are different. We have 

P = Pn + P P 

A Pn - Pp (N-Z) 1 2Z 
Pn + Pp N + Z A 

P P = f(l-A) 

In the Hartree-Fock model at zero temperature, the kinetic energy per nucleon is given 

by 

K.E. _ N3p}{n) Z3p){p) 

A A 5 2m A5 2m 
_ p n 3h / 3p ra . 2 / 3 p p 3h / 3pp , 2 / 3 , . 

p 10m l 8vr ' p 10m l 8vr ' 1 ' 

Expanding the above in powers of A upto A 2 we get 

^ = J_/ l 2(^_)2/3 (1 + 5 A2) ( } 

A 10m v 16vr ; v 9 ; y 1 

This then identifies C s (k.e): 

Of course, in terms of higher powers of A exist which will be small and are neglected. 

For contribution C s (p.e) we start with the simplest potential energy density that will 
produce the correct saturation density, binding energy, compressibility and symmetry energy 
coefficient js}: 

, . A u A[ , 9 o\ B p f v 

V{Pn, Pp) = PnPp +^—{Pn + Pp) + , ] a ( 8 ) 

Po ■tpo o + J- Po 
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This will give for potential energy per particle: 

A ppo 2pp v o + 1 po 

Writing p n , p p in terms of p and A we get 

^ = I > (4, + At) + I > (4 - A„)A* + -£-(.£.)• (10) 
A 4 po 4 po cr + 1 po 

We identify C s (p.e) = \-^{A t - A u ). With A u = -379.2MeV, A l = -334.4MeV, 5=303.9 

MeV and a— 7/6, one gets for symmetric nuclear matter saturation density p = 0.16fm -3 , 

B.E./A=16 MeV, compressibility=210 MeV and at p/p =l, C a (p.e)=11.2 MeV. Together 

with C a (fc.e) « 12.3 Mev C s adds to the total value of 23.5 MeV. 

The Hartree-Fock energy of an orbital is given by 

e = p 2 /2m + A u (p u /p ) + A t ( Pl / p ) + B(p/p ) a (11) 



' 3 ./l! C\P i{( „ - U n )} + 1 



The value of p„ is found by solving for a given p n and (3 = 1/T 

87r r 00 p 2 rfp 
/i 3 7o exp[/3(e n - p n )] + 1 

In the model pursued here, the potential part in e is constant for given densities and is given 

by K = £-{A u + Ai) + ^4g^A + B(f o Y = £(4, + At) + 2C s (p.e)A + For zero 

p 2 

temperature, the chemical potential (at zero temperature p = e/) is p = + X. We first 
consider zero temperature. This will be followed by the finite temperature case. 

The change in neutron chemical potentials for two nuclei: one with A 2 = 1 — ^ and 
another with A 1 = 1 — ^ is 

Sp n = p„(2) - p„(l) = ^M(2) - ^M(l) + 4C s (p.e)(| - |) (13) 
The kinetic term is: 

f?< 2 > - fr' 1 ' - £(|' 2/3 Kf ) 2/3 p + A >> 2/3 - <f ) 2/3 )' 1 + A '' 2/3 i < 14 > 

Expanding the above to the lowest order in A 



2m V ; 2m V ; 6m V 167r % 



(15) 



Together then we get 

« 4 n * 

A 2 



5p n ^A.0*C s (^-^) (16) 



the approximate nature of the above equation arises because contributions from kinetic 
energy have been retained to lowest orders in A. Note the difference of the above equation 
from the generally used relation of eq.(4). 

Consider now finite temperature mean field theory. The contribution to 5fi n from poten- 
tial energy does not change. But the contribution to <5/i n from kinetic energy will change. 
An approximate answer, quite accurate upto 6 MeV temperature, is [-^-(2) — ^-(1)](1 + 

2<"T">2 (^) (ti) 

i2e^(i)e (2) ) w ith eo(l) = 2^~(1) and eo(2) = -|^-(2). But it is easy to get an accurate an- 
swer for all temperatures numerically. Fig.l shows that in this model eq.(16) works better 
than eq.(4). 

One can do refinements to this model. For example as the temperature increases, the 

n q 

nucleus will expand [6|, 17J which will cause some quantitative changes. But we will not 
pursue these finer details. 



IV. 6fj, n IN THERMODYNAMIC MULTIFRAGMENTATION MODELS: CANON- 
ICAL AND GRANDCANONICAL 

We now go back to the multifragmentation model that we briefly alluded to in section II. 
Assume that the system with Aq nucleons and Z$ protons has temperature T, has expanded 
to a higher than normal volume and the partitioning into different composites can be calcu- 
lated according to equilibrium statistical mechanics. In a canonical model, the partitioning 
is done such that all partitions have the correct A , Z Q (equivalently N , Z ). Details of the 
implementation of the canonical model can be found elsewhere jjj; here we give the essentials 
necessary to follow the present work. 

The canonical partition function is given by 

m,j 

fec^En^ (17) 

Here the sum is over all possible channels of break-up (the number of such channels is 
enormous) which satisfy Nq = 1 x n i,j an d Z = J x n i,j'i is the partition function 
of one composite with neutron number / and proton number J respectively and n^j is the 
number of this composite in the given channel. The one-body partition function uji^j is a 
product of two parts: one arising from the translational motion of the composite and another 
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from the intrinsic partition function of the composite: 

WlJ = ^(27rmT) 3/2 A 3/2 x z IyJ {int) (18) 

Here A = I + J is the mass number of the composite and Vf is the volume available for 
translational motion; Vf will be less than V, the volume to which the system has expanded 
at break up. We use Vf = V — V , where V is the normal volume of nucleus with Z protons 
and No neutrons. In this calculation we have used a fairly typical value V = 6Vo- 

The probability of a given channel P(nj t j) = P(n 0i i, ni,o, ^1,1 ni t j ) is given by 

1 / , ni ' J 

P(ni,j) = ^U^ L J (19) 

The average number of composites with / neutrons and J protons is seen easily from the 
above equation to be 

(ni,j) = ^,j Qn ^- J (20) 

The constraints iV = J] / x rij^j and Z = J] J x rijj can be used to obtain different looking 
but equivalent recursion relations for partition functions. For example 

Qn ,z = -rr X] Iui,jQn -i,Zo-J (21) 
7V o 1, j 

These recursion relations allow one to calculate Qn ,z 

We list now the properties of the composites used in this work. The proton and the 
neutron are fundamental building blocks thus zio(int) = Zo^iint) = 2 where 2 takes care 
of the spin degeneracy. For deuteron, triton, 3 He and 4 He we use zi t j{int) = (2si t j + 
1) exp(— f3E It j(gr)) where j3 = 1/T, Ejj(gr) is the ground state energy of the composite 
and (2si y j + 1) is the experimental spin degeneracy of the ground state. Excited states for 
these very low mass nuclei are not included. For mass number A = 5 and greater we use 
the liquid-drop formula. For nuclei in isolation, this reads (A — I + J) 

Zl ,j(int) = exp ±[W A - a(T)A 2 / 3 - - C s ^^- + ^] (22) 

The derivation of this equation is given in several places so we will not repeat the 

arguments here. The expression includes the volume energy, the temperature dependent 
surface energy, the Coulomb energy and the symmetry energy. The term 2^4 represents 
contribution from excited states since the composites are at a non-zero temperature. 



For most of the calculations here the dissociating system is N = 93, Z = 75, A = 168 
for reaction 1. For reaction 2, the dissociating system is iVo = 111,-Zo = 75, Aq = 186. 
These will represent 112 Sn+ 112 Sn central collisions and 124 Sn+ 124 Sn central collisions after 
pre-equilibrium particles are emitted. These two systems have received much attention in 
the past. We also have to state which nuclei are included in computing Qn ,z ( ec i-(17)). F° r 
/, J, (the neutron and the proton number) we include a ridge along the line of stability. The 
liquid-drop formula above also gives neutron and proton drip lines and the results shown 
here include all nuclei within the boundaries. 

The long range Coulomb interaction between different composites can be included in an 



approximation ca^ 
scheme set up in 



led the Wigner-Seitz approximation. We incorporate this following the 

Computations of observables with the canonical model can be done without an explicit 
use of a chemical potential. We can, however, compute the chemical potential using the 



thermodynamic identity /i = {dF / dn)v,T jlQfl- We know the values of Qn ,z ,Qn -i,z and 



Qn ,z -i- Since free energy F is just — TlnQ we compute fi n from — T(hxQN 0t z -^QN -i,z ) 
and /ip from -T(lnQ NotZo -hyQ No>Zo -i). 

We now briefly review the grand canonical model. For finite systems such as considered 
here it is inferior to the canonical model but is easier to implement. If the numbers of 
neutrons and protons in the dissociating system are N and Z respectively, the ensemble 
contains not only these but many others but the average value can be constrained to be 
Nq and Zq. The chemical potentials fi n and fi p serve to fix the average numbers. If the 
neutron chemical potential is fi n and the proton chemical potential is fi p , then statistical 
equilibrium implies that the chemical potential of a composite with N neutrons and Z 
protons is fi n N + fi p Z. The following are the relevant equations for us. The average number 
of composites with N neutrons and Z is 

(n N , z ) = e^ N+ ^ z u N>z (23) 

There are two equations which determine /i n and fi p . 

N = Y / Ne^ N+ ^ z u N , z (24) 



Ze^ N+l3 ^ z u N , z (25) 
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The sum here is over all nuclei within drip lines whose (N, Z) do not exceed (No, Zq) since 
there can not be a composite whose N, Z exceed those of the system from which it emerges. 

We want to point out the following feature of the grand canonical model. In all cu^z's 
in the sum in the above two equations, there is one common value for Vf (see eq.(18)). We 
really solve for N /Vf and Z /Vf. The values of \i n or n p will not change if we, say, double 
No, Zq and Vf simultaneously provided the number of terms in the sum is unaltered. We 
then might as well say that when we are solving the grand canonical equation we are really 
solving for an infinite system (because we know that fluctuations will become unimportant) 
but this infinite system can break up into only certain kinds of species as are included in the 
above two equations. Which composites are included in the sum is an important physical 
ingredient in the model but intensive quantities like [3, \i depend not on N , Z but on N /Vf 
and Zo/Vf. To apply the grand canonical model to finite systems after solving for /x's we 
plug in the value of Vf that would be appropriate for the system N , Z . If the system which 
we are investigating is small, experimental data may show substantial deviations from the 
grand canonical model as we will verify later. 

For later application, we will also use a slightly different version of the above equations 



We label two other chemical potentials: fi (fixes baryon number) and v (fixes total 



charge) : 

A = Y.Ae^ A+ ^ z u A ,z (26) 



Z = £ Ze^ A+ ^ z Co A , z (27) 

Here /i = v = n p — /i and uja,z = &n,z with A = N + Z. 

For A = 168, Z = 75 and A = 186, Z = 75 the /i n 's of the grand canonical model are 
compared with those derived from the canonical model in Fig.2. The values are very close. 
We also compare the 5[i n in the two models. 



V. COMPUTATION OF 5fj, n AST^O 



The computation of whether in canonical or grand canonical requires solving compli- 
cated equations. However it may simplify as T — > (see also [llj]). We try this in the 
canonical model first. 
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As T — > the translational degree of freedom can be considered frozen. Let A , Z be 
stable against spontaneous dissociation (if A , Z is one of the nuclei within neutron and 
proton drip lines then it can not spontaneusly decay into a neutron (proton) plus a daughter; 
usually the only other channel one needs to check is an alpha plus daughter). As T — > 0, 
the system will drop to the ground state of A Q , Z and we will have (eq.(17)) Qn ,z ~^ 
un ,z with uj No ,z now given by exp(-±E gr (N , Z )). Consequently we get fi n (N , Z ) = 
E gr (N , Zq) — E gr (N Q — 1, Zq). This result is of course physically meaningful: —fi n is simply 
the separation energy required to free a neutron from bound nucleus N , Z . Thus 

H n (N , Z Q ) = -W + terml + term! (28) 
terml = a[A 2 /3 - (A - l) 2 / 3 ] + kZ 2 [A 1/3 - (A - I)" 1 / 3 ] (29) 



\N -Z ) 2 (N -l-Z 



For us term2 is of greater significance. It can be rewritten as 

AZ 2 

terml = Cj 1 



\2 



^o^o-l) 7 
A7 2 

« C„(l - ^) (31) 

For 5/i n we need to take the difference between two chemical potentials. Terml of eq.(28) 
contributes little in the difference and thus we end up with familiar eq.(4), i.e., 5/i n = 

4 * Kfr) 2 

Eq.(16) can also be obtained but as an approximation to eq.(30). We can rewrite eq.(30) 

as 

term2 = 2C S - c}-^ - C a {No ~^~^ (32) 

Aq Aq[Aq — I) 

For the examples we are using, A = 168(186) and Z = 75, the third term in the right hand 
side of the above equation is much less important than the second term. If we neglect the 
third term and as before also terml we end up with eq.(16), i.e., 5/j, n = 4* (Zi/A\ — Z2/A2). 

Let us see if we can get a sensible answer in the grand canonical model. We will use 
the the alternative forms eqs (26) and (27): /i controls the baryon number and v the total 
charge. So long as we maintain the general form of eqs. (26) and (27), that is, include in the 
sum all particle stable nuclei with A < A , Z < Z the limits of /i and v are very difficult to 
obtain even at the zero temperature limit. In particular we can not have mutiplicity 1 in the 
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ground state of the dissociating system A , Z Q and zero occupation in all other composites 

n 

in the sum in eqs. (26) and (27) (see, however, [llfl). Exclusive occupation in the ground 
state of A , Z can be achieved in the canonical model but not in the grand canonical 
ensemble as this violates a fluctuation equation. We will deal with that equation but let 
us look at this first in a more pedestrian fashion. At zero temperature the translational 
degree can be considered frozen. From eqs. (26) and (27), to get mutiplicity 1 we need 
(fiA + vZq — E gr (A Q , Z )) — > as (3 — > oo. This alone is not enough to determine /x or v 
but we also require fiA + vZ — E gr (A , Z) to maximise at Z = Z so that at other values 
of Z the difference is negative and occupation in Z's other than Z will go to as (3 — > oo. 



The maximisation condition gives 11 1 



v = -AC S + (33) 

A A ' 

Having determined u, the value of \i can be found from fiA + uZ — E gr (A , Z ) = 0. This 
precudure ensures that fiA + vZ — E gr (A , Z) is negative for Z < Z and hence as (3 — > oo 
the occupation in composites labelled by A , Z with Z < Z will go to zero. But this does 
not guarantee that fiA + vZ — E gr (A,Z) will be less than zero for all A's less than A 
with Z's less than Zq that are in the sum of eqs. (26) and (27). In fact they are not all 
negative and whenever they are postive, multiplicities for those (A, Z)'s blow up. Another 
way of understanding this is to realise that for sole occupation in (Aq, Zq) there are three 
conditions to be met: fiA + vZ — E gr (A, Z) must (1) go to at A , Z , (2) must maximise 
as a function of Z at A , Z and (3) must maximise as a function of A at A Q , Z . With only 
two parameters /x and v this can not be achieved. 

We can also deduce the impossibility of exclusive occupation in the ground state of A , Z 
from very general arguments about fluctuations. It is easy to prove this when there is only 
one kind of particle (eqs. (19) to (21) in |9|]). With 2 kinds of particles, neutrons and protons 
and hence 2 chemical potentials /i and v the notation gets complicated. Quite generally, the 
equation for the grand canonical partition function, when there are many species % which 
are non-interacting, is given by 

oo 

^.can = n(E( e ^T^«) (34) 

i ni=0 

In our case i stands for both a and z, the composite mass and charge; z ni (i) is the canonical 
partition function of rij particles of type i; rii goes from to oo as we are constructing a 
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grand partition function. We have //j = fia + vz. We need not specify the functional form 
of z ni (i). 

Eq.(34) can be rewritten as 

oo 

Zgr.can ^ ' 6 ^ Zm (35) 
m=0 

where we have absorbed the factors e l3uz in z m which is quite complicated but it only contains 
partition functions with total particle number m. We clearly have (fluctuation equation): 

(m 2 )-(m) 2 = ^-^\nZ gr . (36) 

In our case In Z gr is particularly simple: 

\nZ gr = Y,e^ A+ ^ Z Cu A , z (37) 

and 

~^ = J2^M (38) 

If (ua) = 1 for A = A and for all others then (m 2 ) — (m) 2 = A\ but this contradicts the 
assumption that the only occupation is for m = A in which case the fluctuations would 
have been 0. 

In spite of this conceptual difficulty, the expression for \x n derived from the grand canonical 
ensemble in and that derived here from canonical ensemble are not that different. In 
particular 5fi n will be practically the same. For completeness we write the two, one after 



the other. That in 



ll| is 



»n = -W + -^-.J, + C.[l - (-^) 2 ] (39) 



whereas we get 



H n = -W + a[A 2 /3 - (A - l) 2 / 3 ] + kZ 2 [A 1/3 - (A - I)- 1 / 3 ] + C s [l - (^) 2 ] (40) 



VI. 5fi n IN MORE GENERAL CASE 

The contribution of the symmetry energy to Sfi n in the more general case will be function 
of temperature and not a constant as implied in eq.(4). This can most readily be seen by 
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analytically deriving 5fi n at high temperature. Let us derive this in the canonical ensemble 
first. At very high temperature we will get only neutrons and protons and so eq.(17) becomes 
particularly simple: 

, ,N , Z 

Qn ,z - (41) 

The formula \i n = -T[ln Qn ,z ~ ^Qn -i,z ] leads to /i n (N ,Z ) = -Tin ^ and hence 
S/j, n = Tln(N 2 /Ni) and thus not a function of the symmetry energy at all. It is easy to verify 
that in this high temperature limit the grand canonical ensemble gives identical answers. 
From T = towards a large value of T this must happen gradually and so Sfi n must be an 
evolving function of T. 

We are unable to derive a simple formula for 5/i n for a general T. The reasons for 
this failure are obvious enough. Eqs.(24) and (25) are highly non-linear and thus the grand 
canonical ensemble provides no simple leads. The relevant equations for the canonical model 
are easy to compute on a computer but otherwise are quite non-transparent. In Fig. 3 we 
show as a function of temperature (relates to beam energy of collision) numerical results 
for 5fi n for the two systems studied here. Two results are shown; one with C s = 23.5 
MeV and another with C s = 15 MeV. Extracting the value of C s from S/i n using eq.(4) is 
difficult because of temperature and hence beam energy dependence. Notice that at high 
temperature, the curves obtained with different values of C s approach one another. They 
will converge to one line which will rise linearly with T. 

A very interesting aspect of the experimental data is the behavior of 5fi n when Z 2 /A 2 
varies for fixed Z\jA\ or vice versa. In the thermodynamic model this depends on the beam 
energy. There is no such dependence in eq.(4). The simple prediction of eq.(4) is compared 
with thermodynamic model predictions in Fig.4. In the model, for fixed Zi/Ai, 5/i n is still 
approximately linear with (Z 2 /A 2 ) 2 although deviation can be seen by eye-estimation. It is 
also noticed that use of eq.(4) to estimate C s from 5/i n will underestimate the value of C s 
at low temperature (4 or 5 MeV) but will overestimate the value at high temperature, for 
example, at 7 MeV. 

Experiments directly measure a — (3fj, n rather than /i n and we plot a from the ther- 
modynamic model as a function of temperature in Fig. 5. We also obtain a from eq.(4) 
and compare. Eq.(4) gives a 1/T dependence but the fall with T is much slower in the 
thermodynamic model at higher temperature. In the thermodynamic model a would reach 
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asymptotically a constant value ln-j^-. This difference in behavior between the two predic 
tions can be ascertained in experiments. 

There has been much activity in recent times relating a to C s Q, Q, Q] . An approximate 
derivation of eq.(4) from the expanding excited source (EES) model can be found in [3 



Attempts to obtain this from antisymmetrised molecular dynamics can be 
Temperature dependence of symmetry energy was discussed by Li and Chen 



bund in |l5] 
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VII. DEVIATION FROM ISOSCALING: EXPLANATION FROM CANONICAL 
MODEL 

We now discuss isoscaling and possible deviations from it. It will become clear in this and 
the following section that, from the point of view of theory, isoscaling can work very well 



when N/Nq and Z/Zq are small (< 0.35). Many experimental data fall in this range 
and isoscaling is one robust feature to emerge in recent experimental intermediaite energy 
heavy ion collisions. If we now extend these observations to larger composites deviations are 
to be expected. Will the grand canonical and canonical results for i? 21 always agree? The 
answer is no; it depends upon the size of the dissociating system (N ,Z ), more precisely, 
upon the fractions N/N Q and Z/Z . The grand canonical model always predicts isoscaling. 
In this model 

p _ „(r^z(2)) 

-«21 — 7TTT 

(njv,z(l)> 

= Cexp(aN + pZ) (42) 

where we have used eq.(2) and the advantage that with similar beam conditions the factors 
cl>tv,z in the denominator and the numerator cancel each other out. Therefore, in the grand 
canonical approximation, the slopes of In R21 as a function of N for fixed Z (a) will never 
deviate from a straight line and (b) for different fixed Z's the slopes will not change. In 
many experiments where the sizes of the composites fl encompass from small to large this 
is not true. In Fig. 6 we compare some experimental data with a grand canonical calculation. 



Experimental details can be found in [J, [17, 18]. In the experiment, reaction 1 is 58 Ni on 
9 Be and reaction 2 is 64 Ni on 9 Be. For the grand canonical calculation, for reaction 1 the 
dissociating system is taken to be 58 Ni+ 9 Be (AT = 35, Zq = 32) and for reaction 2 the 
dissociating system is taken to be 64 Ni+ 9 Be (No = 41, Zq = 32). All composites between 
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drip lines are included as detailed in section IV with the highest values of N, Z terminating 
at N , Z . The limitations of the grand canonical model are very obvious in the figure. The 
slopes of In R 2 \ differ from straight lines for large Z and as well the same slope for all Z 
does not fit. In these calculations the same temperature was assumed for all composites. 
The actual situation may be more complicated requiring a different temperature for higher 
composites. However the deviation from linearity will require even further complications 
if one insists on using the grand canonical model to fit the data. From the point of view 
of theory, however, for the emitting systems in these cases, the use of the grand canonical 
approximation for the emissions of heavier composites is not valid. 

The canonical model does not impose these restrictions. Now (see eq.(20)) 



This formula is not transparent at all but produces deviations from isoscaling for small 
systems. 

Let us show the results of the canonical model calculation for the same case as above: Ni 
on Be. The parameters for the calculations are the same as used for the grand canonocal 
model: the same temperature, the same freeze-out volume and the same composites included 
in building the respective partition functions. But the canonical calculations (Fig. 7) are sig- 
nificantly different from the grand canonical results (Fig. 6) and much closer to experimental 
data. 



VIII. FURTHER OBSERVATIONS ABOUT DIFFERENCES BETWEEN THE 
CANONICAL AND GRAND CANONICAL RESULTS 



relating (n N>z ) in canonical with that in grand canonical. The formula for (n^z), ec f(20), 
respects constancy of total neutron and proton numbers in the partition function and in this 
regard is superior to eq.(23). If N/N or/and Z/Z are not small (for some cases in Figs. 6 
and 7 they are > 0.5), eq.(20) can give quite different results from those given by eq.(23). 
This is so even when the chemical potentials given by the two models are very close. If the 
chemical potentials in the two ensembles are nearly the same (as found in all our examples 




(43) 



Instead dealing directly with i? 2 i oc 




' z z \i^ let us investigate a simpler task: that of 
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so far) that merely guarantees that 



Qn -i,Zq _ ^ 
Qn ,z 



Qn ,z 

Here the left hand side is given by the canonical model (and leads to the definition of fugacity 
and chemical potential in the canonical model) and the right hand side is computed from 
the grand canonical model (eqs.(24) and (25)). However, for {n^,z) to be nearly the same 
in the two ensembles we require 

Qn -n,Zq-z ^ e PnnN+PnpZ £41^ 
Qn ,z 

We show below how the left hand side of the above equation can lead to the right hand side 
when N/No and Z/Zq are small but are expected to deviate when they are not small. To 
proceed let us call ^q' 1 ^ =neutron fugacity of the system A, B and Q q^ b 1 =proton fugacity 
of the system A, B. Rewrite the left hand side of eq.(45) as 

Qn -n,z -z _ Qnq-n,z -z x Qn ,Zq-Z 
Qn ,Z Qn ,Z -Z Qn ,Z 

The first term on the right hand side can be rewritten as a product of N terms involving 
neutron fugacities and the second term as a product of Z terms involving proton fugacities. 

Qn -N,Zq-Z _ Qn -N,Zq-Z Qn -N+1,Zq-Z Qnq-1,Z -Z 
Qn ,Z -Z Qn -N+1,Z -Z Qn -N+2,Z -Z Qn ,Z -Z 

The first term in the right hand side of the above equation is the neutron fugacity of 
the system N — N + 1, Z Q — Z, the second term is the neutron fugacity of the system 
Nq — N + 2, Zq — Z and so on; finally ending with neutron fugacity of the system No, Zq — Z. 
If TV is negligibly small compared to iVo and also Z is negligibly small compared to Zq then 
each of these terms can be approximated by neutron fugacity of the system N , Z (for iV w 
20 and Z w 20 in Fig. 7, this leads to gross errors) leading to ( Q ^°' 1 ' z ° ) JV . Equating %a^Mfi 

^N ,Z ' 1 ° Qn ,Z q 

to e^™ we get the factor e^™^ of eq. (45). (Even this can introduce significant error for 
iV w20.) It is clear also how e^ z can arise by resolving Q ^°' z °~ z into Z proton fugacities. 
The connection between the two sides in eq. (45) can also be established using a saddle- 



point approximation 19| but we will not pursue this any further. 

It then follows that although at low N, Z canonical and grand canonical calculations can 
agree they will diverge when N, Z grow. This is highlighted in Fig.8. 
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IX. DISCUSSION ABOUT SECONDARY DECAYS 



Before comparing with data (as is done in Fig. 7) one needs to investigate the effects of 
sequential decay on calculated i?2i- The multiplicities (n^j) refer to populations of nuclei at 
finite temperatures. Nuclei at finite temperatures can decay by particle emissions and values 
of {ni t j) will change. However, because one is comparing ratios, the effect on R21 may be less 
drastic. We have included the contributions from excited states through a factor T 2 A/eo in 
eq.(22). This overestimates contributions to the partition function from excited states and a 



cut-off will be necessary 



20j. This requires detailed work which we have not carried out. In 



investigating isoscaling using antisymmetrised molecular dynamics, Ono et al. 21[ find that 
the effect of secondary decay is to decrease a to a/2. This means the experimental value 



of logR,2i should be compared to logyn(2)/n(l) rather than to log[n(2) /n(l)]. We have 
applied this "empirical" correction to the Ni on Be case in Fig. 9. We now use temperature 
T=5 MeV. This leads to a steeper rise (compared to the T=8 MeV). But after the correction 
the rise decreases to a value more compatible with experiments. It now corresponds roughly 
to "uncorrected" T=8 MeV calculation. 



X. SUMMARY 

This work has addressed two issues. How well does eq.(4), a much used formula, fit results 
from calculations obtained from well-known models of nuclear dissociation? We find the fit 
approximately correct though not exact (the mean field model gives a different formula but 
the mean field model is not a model for multifragmentation). A more important question 
is: how do we handle cases when deviations from isoscaling are significant. Preliminary 
results reported here are very encouraging. Statistical multifragmentation can explain the 
results but one must give up using the grand canonical ensemble and instead use models 
that strictly conserve neutron and proton numbers. 
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FIG. 1: Mean-field model calculation (dashed curve) for 5/U ra compared with simple versions 
(eq.(16) and eq.(4)). 
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FIG. 2: Chemical potentials fi n (top panel) for systems Z\ = 75, Ai = 168 and Z2 = 75, A2 = 186. 
Canonical model (dashed curve) and grand canonical model (solid line) results are quite close. The 
bottom panel shows results for 5fi n in a much more enlarged scale. 
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FIG. 3: The difference Sfi n calculated from theory compared with the simple prediction of eq.(4). 
Since canonical and grand canonical results for 8/j, n are so close we just display the canonical model 
results. Two values of the symmetry energy coefficients were used. Notice that thermodynamic 
model predictions for 5fj, n depend on the temperature unlike the simple model prediction (eq.(4)). 
In particular, irrespective of C s , the theoretical value of 5fi n will approach the value T In N2/N1 . 
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FIG. 4: This figure tests the linearity of <5/x n against (Z2/A2) 2 for fixed (Zi/A±). Canonical model 
calculations are done for varying values of Z2/A2. The calculations are done with the value of 
C s at 23.5 MeV. The solid line is a plot of eq.(4) with (7 S =23.5 MeV. Notice that at 4 or 5 MeV 
interpreting the canonical model result using eq.(4) will lead to un underestimation of the value 
of the underlying C s and an overestimation at 7 Mev temperature. Linearity in thermodynamic 
calculation is not perfect but roughly correct in the temperature range displayed. 
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FIG. 5: The isoscaling parameter a = 5fi n /T calculated in the canonical model (the dotted curve) 
is campared with the standard version 4(C S /T)[(Zi / A\) 2 — (Z2/A2) 2 ] (solid curve). The thermody- 
namic model predicts that a becomes constant asymptotically whereas in standard parametrisation 
this would fall off like 1/T. 
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FIG. 6: Ratio of cross-sections of producing the nucleus (N,Z) where reaction 1 is 58 Ni on 9 Be and 
reaction 2 is 64 Ni on 9 Be, both at 140 MeV/n beam energy. Experimental data with error bars are 
compared with theoretical results from grand canonical ensemble ( hollow points). Dotted lines 
are drawn through experimental points and solid lines through calculated points. For calculation 
the dissociating systems are taken to be Ni+Be. The constant C (eq.(42)) should be of the order 
1 and in drawing the figure is taken to be 1. The value of the constant does not affect the slope of 
the log of the ratio. Notice that for large (N, Z) the slopes are no longer those of a straight line 
and that "average" slopes for small (TV, Z) and large (N, Z) are very different in experiments but 
the same in grand canonical calculation. 
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FIG. 7: The same as in Fig.6 except that the calculation is done in the canonical ensmble. Agree- 
ment with data is far superior compared to that in Fig.6. 
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FIG. 8: Comparisons of R21 in canonical and grand canonical models. The solid lines are through 
points computed in the canonical model and the dotted lines are through points calculated in the 
grand canonical model. The dissociating systems are the same as in Figs. 6 and 7 and most of the 
points also appear in those figures. Here we show results for lower N, Z also to emphasize that for 
low N, Z agreement between canonical and grand canonical is very good (and for these isoscaling 
works) but that they differ widely for large N, Z. The grand canonical model should be discarded 
for large N, Z. 
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FIG. 9: The solid lines are canonical model results for R21 at temperature 8 MeV (no effects of 
secondary decay). This already fits the data quite well (Fig. 7). For temperature 4.5 MeV we show 
two results. The dashed lines are the direct calculations (no corrections for decay). The slopes 
are much higher than the 8 MeV curves. The dotted curves include a phenomenological correction 
(see section IX) due to decay. It brings down the slope and mimics the raw calculation at 8 MeV 
(solid line here) and hence would fit the data reasonably well. 
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